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By an extension of the Bethe ansatz method used by Gwa and Spohn, we obtain an exact expression 
for the large deviation function of the time averaged current for the fully asymmetric exclusion 
process in a ring containing N sites and p particles. Using this expression we easily recover the exact 
diffusion constant obtained earlier and calculate as well some higher cumulants. The distribution 
of the deviation y of the average current is, in the limit N — > oo, skew and decays like exp — (Ay 5 ' 2 ) 



OS , for y — » +oo and exp — (A'|y| 3//2 ) for y — > — oo. Surprisingly, the large deviation function has an 
OS ■ 

, expression very similar to the pressure (as a function of the density) of an ideal Bose or Fermi gas 

|. in 3d. 

jy^ ! PACS:02.50.+s, 05.40.+j, 82.20-w 

In the last decade, it has become clear that questions as different as Id turbulence, growth of interfaces or directed 
t-H , polymers in a random medium in dimension 1 + 1 dimension were different versions of the same problem described by 
t^J- ■ the noisy Burgers or equivalently the Kardar-Parisi-Zhang (KPZ) equation 0,^]. One of the simplest lattice versions 
of this problem is the asymmetric exclusion process (ASEP) which also has a long history both in the mathematical 
and in the physical literature [^^|. It describes a driven lattice gas with hard core exclusion and under suitable 
scalings the evolution of its macroscopic density is governed by a noisy Burgers equation in Id which is equivalent to 

OS . 

the KPZ equation. 

Many properties of the ASEP have been calculated (steady states, mass gaps, diffusion constants) by various 
methods including exact solutions |]6|-|l3|| and compared @JlJ] with the predictions of the replica approach |l5|-|l7|l or 
Ch , numerical simulations |fl8|-pof for directed polymers or growth problems. 

Here, extending the Bethe ansatz approach of Gwa and Spohn [p|, we obtain the expression for a large deviation 
function which gives the exact probability for the time averaged current to take any, typical or untypical, value. 
Similar distributions have been studied numerically in the past []l8|-pc| , in particular to better understand the n — > 
limit in the replica approach of directed polymers in a random medium. Other recent results include the numerical 
calculation of the distribution of the density for a randomly driven Burgers equation plj] . The results presented below 
are, to our knowledge, the first exact determination of a non-trivial probability distribution for a fluctuating quantity 
described by the noisy Burgers equation or equivalently the KPZ equation. Through the known equivalences between 
the different problems mentioned above, our results could easily be used to obtain the whole distribution of the ground 
state energy of a directed polymer in a random medium in dimension 1+1 for a finite lateral geometry. 

We consider here a totally ASEP consisting of p particles on a ring of N sites. During any time interval dt, each 
particle jumps with probability dt to its right if the target site is empty. The unique stationary state of this system 
is one in which all (^) different configurations C have equal weight ||. We call Yt the total distance covered by all 
the particles between time and time t (thus Yq = 0). The probability Pt{C) of finding the system in configuration 



1 



C at time t satisfies 

jP t (C) = £[M (C, C) + Mi(<7, C")]P*(C") 
c 

where Mi(C, C')dt is the probability of going from C to C during a time interval di, increasing Y by 1, while M is 
a diagonal matrix, 

M (C,C) = - J2 M 1 {C',C) 

Let Pt(C,Y) be the joint probability of being, at time t, in the configuration C and having Y t = Y, some non- 
negative integer. Define also F t {C) as ^y =0 exp(aF)P t (C, F). Then 

= £[M (C, C") + e^M^C, C')]F t {C) 

Summing F t (C) over C we obtain (e aYt ). It is easy to see that for large t 

(e aYt ) ~ e x(a)t (1) 

where A(a) is the largest eigenvalue of the matrix Mq + e a M\ (we have omitted the dependence of A on p and N). By 
~ we mean: take the logarithm of both sides of (1), divide by t and let t — > oo to obtain an equality. While (1) holds 
for every initial configuration, it is sometimes convenient to think of the initial state Po(C) as the stationary one. 

In the present paper we show that the exact expression of A(a), for arbitrary p and N, is given by the power series 
in a obtained by eliminating B between the following two series 

^--- P p \ P JZ-7 q -,y. (2) 

^ (pq)l (Nq —pq)\ W 

These series have a non-zero radius of convergence in B so A(a) is analytic in a at least near a = 0. Moreover, by 
the Perron- Frobenius theorem there is no crossing of the largest eigenvalue of the matrix M + e a M 1 for real a and 
fixed finite N. Hence (§) and (j^) determine \(a) on the whole real axis as the analytic continuation of the expression 
in the vicinity of a = 0. In the rest of this work, we first show how (||, |^) can be obtained using the Bethe ansatz. 
Then we investigate some consequences of this exact expression. 

The Bethe ansatz we use is very similar to that used by Gwa and Spohn ||, except that here there is an extra 
parameter a. Specifying a configuration C by the positions 1 < x\ < X2 < ■■■ < x p < N of the p particles we write an 
eigenvector of the matrix Mq + e a M\ in the form 

£^E[[*rf' (4) 

Q 3=1 

where the sum is over all the permutations Q of 1, 2, ...,p and the p parameters Zk satisfy the following equations 

{ea ~ N Zk)P =(-) p - 1 f[(e a -z j ) (5) 
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A(a) = e a 



P (7) 



z 1 z 2 .-z p = 1 (6) 

The derivation of (^||) is almost identical to that given by what Gwa and Spohn in Q and we refer to their work. 

For any solution {zu} of (||J1), expression dJ) gives an eigenvector of the matrix M + e a M\ with eigenvalue X(a) 
given by 

"11 1 " 

1 h ... H 

zi z 2 z p _ 

The difficulty is that there are in general many solutions of (^^) corresponding to different eigenvalues of the matrix 
AIo + e a M\. The one we need to choose is the largest eigenvalue A(a) which, as already noted, can be followed when 
a varies, from the unique eigenvalue which vanishes when a — > 0. Thus A(a) corresponds to the solution such that 
all the Zk — > 1 as a — > 0. 

To see why (|^,^) are equivalent to (|5]|],[|]) when one chooses the above solution, we consider a function h(z) analytic 
in the circle \z — e a \ < 2e. It is then an easy consequence of the Cauchy theorem that 

r ,, 7l I dz ^z-e^lpz-Njz-e")} 

where the contour is the circle \z— e a \ = e, the parameter A is such that \ A\ -^i e p e~ Na and the Zk are the p solutions 
of 

(z - e a f = Az N (9) 

which tend to e a as A — > 0. As (||) has a form identical to (jE)), one can use (||) to calculate A(a) (by choosing 
h(z) = —1 + e a /z) and this gives (|^) in terms of B = Ae^ N ~ p ^ a . The constant B can then be determined by requiring 
that (H) is also satisfied. To do that one can use again the identity (||) for the function h(z) = logz and this gives (|^). 
From the exact expressions (|^||), one can expand X(a) in a power series in a. In particular, one can check that 

p(N-p) 

Pip!) 2 (N-p)l 2 (2JV-3)! 2| ^ 3 



and using (||), this gives 



(2p!) (iV-1)! 2 (27V-2p-l)! 

_ p(N - p) 



a z + 0{a A ) (10) 



ton {Yt) - dA(a) 



AT - 1 

a— U 



(ii) 



which is known to be a simple consequence of the fact that in the steady state all configurations have equal weight 
HQ. (Note that if Pq(C) was the stationary measure, (|ll|) would in fact be valid for all t.) The expansion ( ^fj| ) also 
gives the diffusion constant 

(y t 2 ) - (y t ) 2 d 2 A(«) 



lim 



(12) 

a=0 



t^oo t da 2 

which agrees (up to a trivial factor p 2 due to a slight change in the definition of Y t ) with the expression obtained by 
an extension |Tcj of the matrix method || . Exact expressions of higher cumulants of Y t may be obtained in a similar 
way from (p]|3]q) by expanding A(a) to higher orders in a. For example, the third derivative of A at a = gives 
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Um <y3)-3(iT(y) + 2(y)3 = (13) 

t^oo t 

3p(N-p) 3 (2N - l\ 2 (N- 1 



(2N- 1)(N- 1) \ 2p J \ p 

4p(N~p) 3 ( 3N — 1\ / N — 1\~ 3 
"(3JV-l)(JV-l)l, 3p A p J 

For large positive a (i.e. a>l), the solutions zu of (H]§), for 1 < k < p are given by 

z k = e^-P-W +0(e- a ) 

and this leads to 

For large negative a (i.e. —a ^S> 1), one of the Zk becomes large (say Z\ ~ e^ 1 ^") and the p — 1 other 2* become 
small (zfe ~ e a for 2 < fc < p). This leads to 

-1 + eP" for p<N-p 
X(a) ~ ^ -l + 2eP Q for p = N-p (15) 

-l + e (JV-p)a for p> N- p 

(We can understand the limit a — > — oo in ( |l5| ) by noting that the probability of not making any jump during the 
time interval t is, for a configuration in which all particles form one cluster, just exp(— t).) 

In the limit of large N and p, keeping the ratio p/N = p fixed, one finds for the first four cumulants of Y (by 
calculating the first four derivatives of X(a) at a = 0) 

(Y) 

lim V " N P0- - P) ( 16 ) 



lim (y2) 7 (F)2 * iV^[p(l - p)]>/^ (17) 



(y 3 ) -3(y 2 )(r) + 2(F) 3 

lim 

w_ 

2 3 3 / 2 



N 3 [p{l- P )f[l-4j-\, (18) 



,. (r 4 ) - 3(r 2 ) 2 - 4(r 3 )(r) + i2(y 2 )(r) 2 - q{y) 1 

hm 



N^W~p)? ,2 ( 1 4 + 4b-^)- V2 (19) 



t^oo £ 

2 2 L - :V- - 

In fact, in this limit N — ► oo with p = p/iV fixed, the expressions (|^) and (^|) can be rewritten as 

p(N - p) 



I 2nN 3 
P(l - P) 



A (a) — a 



iV 



- E Hv 2- (20) 
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where 

aV2np(l-p)N^ ~ (3 = £ 1 \ /2 (21) 

q>l q 

and C = -BN N p-P(N - p)- N +P. This shows that for N large and (3 fixed, i.e. for a ~ TV" 3 / 2 , the difference 
A(a) — ap(N — p)/N takes a scaling form 



(22) 

The two sums which define G(/3) in the right hand side of ( J20|]21| ) are very reminiscent of the pressure and of the 
density of an Ideal Bose or Fermi gas (depending on the sign of C) in 3 dimension p2] | : positive C and a correspond 
to the Fermi gas and negative C and a correspond to the Bose gas. 

The range _/V~ 3 / 2 <C a <C 1 (i.e. large positive C) correponds to the high pressure or low temperature limit of a 
Fermi gas and this analogy leads to 

X(a)-a P -^^^ [p (l-p)]^Na^ (23 ) 

Negative a can also be analysed: when C varies from to —1 in ( p(i| , pl"| ), a decreases (this corresponds to a Bose gas 
in its high temperature phase) and a has a finite limit ao as C — ► — 1. To decrease further a, it turns out that one has 
to increase C from — 1 to and to replace the expressions ( pf]| ) and ( pl| ) of G{j3) and (3 by their analytic continuations 
11 

Gift) = |v^hlog(-C)] 3 / 2 + (24) 



9>1 



= -4v5F[- log(-C)] 1 / 2 + £ ( j 5/2 ° (25) 

9>1 9 

This analytic continuation is a major difference with the usual Bose-Einstein condensation where C remains fixed at 
the value —1 in the whole low temperature phase. Here, for a < ao, one of the Zk becomes much larger than the 
others but, in contrast to the Bose condensate, the chemical potential C in (p3) keeps varying with a. For the range 
TV- 3 / 2 < - a < 1, this leads to 

X(a) ~ a * ±W p)]»N*e? (26) 

The knowlegde of the function A(a) determines (through y — d\/da—p(N —p)/N and f(y) = X — a[y+p(N—p)/N}) 
the large deviation function f(y) defined here as 

f(y) = Hm j log jprob - Np(l - p) = y^j J (27) 

From (p^,p^,p2],p3p(j) one can show that for N large, f(y) takes the scaling form for \y\ <C N 



■kN 3 \p(l-p) 

where the function H has the following properties 



/frW^'f^l (28) 



5 



H(y)~-(y-lf for 



(29) 



H(y) ~ -(2V3/(5^) y 5 / 2 for v -f +oo (30) 

~ -(4v^tt/3) |y| 3/2 for y ^ -oo (31) 
Note that the difference between (Y t )/t — p(N — p)/(N — 1) and iVp(l — p) used in (B7j) is responsible for the fact 



that the maximum of f(y) is at y = 1 as in (|2^). The scaling form ( p8[ ) was obtained ( 20 - 2l|) for a ~ N~ 3 ^ 2 and this 
means (p^-pl]) that it is valid for y of order 1. 

For y of order N, one expects from (|14|Jl5|) the following two scaling forms for f(y): 



f(y)~NH+(yN- 1 ) for y > (32) 

f(y) ~ (yA" 1 ) for y<0 (33) 

As y/(Np) is the deviation (from its most probable value 1 — p) in the velocity of the particles, the extra factor N in 
( |32| ) has a simple interpretation: to speed up the traffic, one needs to speed up all the particles whereas to slow down 
the traffic, it is sufficient to slow down a single particle. Note that if one considers that the scaling forms ( ^ , j30|]3l| ) 
remain valid for 1 <C \y\ -C N, it is possible to obtain from ( [30|]3~j| ) the small z behavior of H+(z) and H-(z). 

The distribution of y is skew ( |T^Jl^ ) and for large y (but still \y\ <§; N), it decays ( |3b"| , j3l| ) as the exponential of 
a power law with an exponent 5/2 for y — > oo and 3/2 for y — ► — oo. Skew distributions with similar tails have 
been observed in numerical simulations fl|,|l8|-p0|| and partly understood by the replica approach A precise 

comparison with existing numerical data cannot however be done because in the present work, we consider a finite 
geometry (a ring of N sites, taking the long time limit first, and then the limit of an infinite lattice (N — > oo)) 
whereas the simulations had been usually done for directed polymers on an infinite lattice jjj (thus taking the limits 
in reverse order). It would certainly be interesting to try to extend the above Bethe ansatz calculation to this case. 
Other extensions could be the calculation of the diffusion constant |24j , the higher cumulants of Y or the whole large 
deviation function in more general cases such as the partially asymmetric exclusion process |24|-|27j or the case of 
several species p8) . 

The similarity of the expression ( pO| , pTl ) of G((3) as a function of (3 with that of the pressure of a Fermi or a Bose gas 
in 3d was a surprise to us and it would be nice to have a more direct explanation for it. A Bose-Einstein condensation 
for the ASEP has already been found in the presence of quenched disorder |2{J. Here, the mechanism seems to be 
different as instead of observing the Bose-Einstein condensation, one finds (|2^ , ^5| ) that G(/3) is analytically continued 
beyond the critical (3q. 
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